Finite E <g> /3 Jahn- Teller Systems: 
a Continued-Fraction Approach 



K. ZiegleiB 



Institut fiir Physik, Universitat Ausgburg, D-86135 Augsburg, 



Germany 



A recursive method is developed to treat electrons coupled to phonons. It is applied to small 
systems with E ® f3 Jahn- Teller coupling. Two cases are considered, a model with one electron and 
two orbitals on a single site (related to the Rabi Hamiltonian) and a model with two electrons on 
two sites. The corresponding Green's functions are represented by rational functions. It is found 
that the spectra change substantially when one phonon couples to the electron but are relatively 
robust under an increasing number of phonons. 



Small systems, like molecules or clusters of atoms, have attracted much attention in recent years 
because of new experimental techniques that provide detailed information of their spectral properties^. 
There are two mechanisms that control the physics of these small systems, one is the tunneling of electrons 
between different orbitals and different atoms, the other is the coupling between electrons and vibrational 
modes (phonons) of the molecules or clusters. For the latter the Jahn- Teller coupling scheme is relevant. 
Small Jahn- Teller systems can also be understood as building blocks for lattice Jahn- Teller systems which 
play a crucial role in solid-state physics, for instance, in the form of transition metal oxidsi 

The main problem of treating electrons that couple to phonons is that even for small systems with 
one or a few electrons the Hilbert space is infinite dimensional. This implies a complex spectrum with 
level crossing and avoided level crossing^. There are various treatments of small Jahn- Teller systems, 
e.g. exact numerical diagonalization with truncated phonon spectrum^ii, Monte-Carlo simulations^ or 
variational methods^. In this paper a systematic recursive procedure for treating the phonons in small 
electronic systems is developed and applied to several examples. It is based on a projection formalism 10 . 
The method is quite flexible, can easily deal with degeneracies and was previously introduced to two- 
component bosons on a latticed. The central idea is to approximate the elements of a Green's function 
systematically by standard (e.g. rational) functions. 

The paper is organized as follows: In Sect. 2 the evolution of a quantum state and its connection to 
a projected Green's function is briefly discussed. The continued-fraction approach is developed in Sect. 
3. Then in Sect. 4 two models with Jahn- Teller coupling are introduced, a single-site model where an 
electron tunnels between two orbitals and a model where two electrons tunnel between two sites. The 
application of the continued-fraction approach to these models is explained in Sect. 5 and the results are 
discussed in Sect. 6. 



where the Hamiltonian H is measured in units of %. A Laplace transformation for a positive time t gives 
for Imz < a resolvent: 



Suppose that the initial state |^o) is from a restricted Hilbert space of low energy, the Green's function 
(z — H)^ 1 acts on restricted Hilbert space, represented by the projector Po and the projected resolvent 
(z — H)~ 1 P$. To evaluate the probability for the system to return to the initial state l^o), the following 
quantity must be evaluated: 
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If Ej are the eigenvalues of H we can write for Eq. 

(*o|(,-i/r 1 |*o)=E K f L ^- 

3 3 

The poles of this expression are the characteristic frequencies of the evolution, starting from the projected 
Hilbert space and returning to it. The imaginary part of the Green's function gives the corresponding 
spectral density. 



III. PROJECTION FORMALISM AND CONTINUED-FRACTION REPRESENTATION 



After this preparation the goal is to evaluate the projected Green's function Pq(z — H) 1 P , where Pq 
projects the states of the entire Hilbert space to the subspace 7io- It satisfies the identity 

Po(z - ff)- x Po = [Pq(z - H)P - P HP x {z - H)^P X HP^ \ (2) 

where Pi = 1 — Pq projects onto the Hilbert space Hi that is complementary to Ho- If H satisfies 

P HP 1 = P a HP 2 , P 1 HP a = P 2 HP (P 2 ^Pi), 

Eq. @ becomes 

P (z - HyiPo = [P (z - H)P - P HP 2 (z - H^ 1 P 2 HP a \\ (3) 

The identity used in Eq. J5J) can be applied again to P 2 (z — H)^[ 1 P 2 on the right-hand side. In general, 
a hierarchy of projectors Pk onto Hilbert spaces Hk is produced recursively. It is based on the fact that 
the projector P 2 j+i is created from P 2 j-i and P 2 j as 

H 2 j+i = T~t 2 j-iXH 2 j c H 2 j-i, 

and P 2 j+ 2 comes from the relation 

P 2j+ iHP 2j = P 2j+2 HP 2j and P 2j HP 2j+1 = P 2j HP 2j+2 . 
This construction implies a recursion relation that connects pairs of projected Greens's functions: 



P 2j (z - n)- 2 l_ x P 2] = [P 2j (z - H) 2j _ x P 2j - P 2j HP 2j+2 (z - H)^ +1 P 2j+2 HP 2j 
By using G 2j = P 2j (z - H) 2j 1 _ 1 P 2 j and H jtJ+1 = P 2j HP 2j+2 this reads 



2;; 



G 



2,1 



z - P 2 jHP 2 j - Hjj +1 G 2j+2 H j+ljj 



(4) 



Further simplifications are possible if it is assumed that the Hamiltonian H can be written as a sum of 
two Hamiltonians as H = Hq + H± with the following properties: 

(1) H must stay inside the projected Hilbert space: HqP = P a H P and Po^o = PqHqPq- 

(2) H\ maps from 7i 2 j to H. 2 j+ 2 : 

Hi ■ H 2 j — > 7Y2j+2, 

where H 2 j is orthogonal to H 2 j+ 2 . Examples shall be discussed subsequently. 
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IV. THE MODELS: E ® (3 JAHN-TELLER COUPLING 

The Jahn- Teller electron-phonon coupling describes electrons that may occupy two possible orbitals. 
The latter can be understood formally as a pseudospin with values a =f , \. On a molecule, a cluster or 
a lattice one phonon mode couples to electrons at each site, where the interaction enery is a sum over all 
sites j = 1, 

H ep h = 9 + b j)( n 3l ~ n 3l)- 

J 

rij a is the electronic number operator at site j and orbital a, and bj (bj) is the creation (annihilation) 
operator of a phonon. 

The simplest case considered in this paper is a single site with two orbitals, where an electron can 
tunnel between the two orbitals. This problem is also well-known in atomic physics under the name of 
Rabi Hamiltoniani^. An extension is a two-site system with two orbitals per site and with two electrons. 
The model is simple if both electrons are in the same orbital because of Pauli blocking. Therefore, we 
consider the case where the two electrons sit in different orbitals. They can tunnel between the two sites 
and are coupled to each other via the electron-phonon interaction. 

A. Single site: Jahn- Teller effect with interorbital tunneling 

There are two electronic states, either the electron is in orbital f or in orbital |. Moreover, there are 
N (N — 0, 1, ...) phonons. Thus the Hilbert space is spanned by 

|JV,T>, \N,i). 

For electronic tunneling rate t between the two orbitals the resulting Rabi Hamiltonian reads 

H = ta\ + uoaawb + ga?,(b^ + b), 
where the Pauli matrices refer to the electronic states. 



B. Two sites: Jahn- Teller effect with intersite tunneling 

The extension of the model to two sites and two electrons provides the opportunity to study the effect 
of electron-electron interaction and electronic correlations. For a given pair of electrons, represented 
by integer numbers n Jcr = 0, 1, states with two electrons that occupy different orbitals are considered. 
(Electrons stay in their orbital because inter-orbital tunneling is not included.) Then the following 
electronic states are available: 

I U>, I U>, I IT,o), |o,|T). (5) 

The tunneling of the electrons and their Coulomb (Hubbard-like) interaction is defined by the Hamiltonian 
H e and by the dispersionless phonons with energy u>q as 

H e = -t ( c L c 2<r + cL c 1t) + U Y n n n jl> H ph = h ] h 3- 

<r=T,i J=l>2 3=1,2 

Cj (cj) is the creation (annihilation) operator of a electron. The electronic spin is not taken into accout 
here, i.e. spin-polarized states are considered. The electron-phonon interaction reads 

H eph =9^ ( & j + M("JT _ 
3=1,2 

leading to the total Hamiltonian H = H e + Hph + H ep h ■ 
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V. APPLICATION OF THE CONTINUED-FRACTION APPROACH TO SMALL SYSTEMS 



If Hi is a small perturbation to Hq the resolvent of H = Hq + Hi can be written as a Neumann series: 
(z - H)- 1 = (z-H - Hi)- 1 = {z- Ho)- 1 Y} H ^ Z - Hoy'}'- 

l>0 

Truncation after a finite number of terms yields poles only from the zeros of z — Ho- This is often 
insufficient to observe a realistic pole structure of the Green's function (z — H) . The approximation can 
be improved by using a Pade approximation^ or a partial summation of infinitely-many contributions^. 
A systematic approach is the continued fraction of Sect. 3 that approximates the projected Green's 
function Pq(z — H)~ 1 Pq by rational functions with a complex pole structure. Depending on the regime 
(weak or strong electron-phonon interaction) there are two different approximation schemes. 



A. Strong Electron-Phonon Interaction 



If the tunneling energy is small in comparison with the electron-phonon interaction it is possible to 
separate the Hamiltonian as H = Ho + Hi with 



Ho = LUo&o^b + ga^Qr + b), H\=tai 



(6) 



and consider Hi as a perturbation in the sense of the discussion in Sect. 3. To diagonalize Ho, the 
Lang-Firsov transformation can be used as a unitary transformation 



e -^-b) } ' 



= .gM) 



with 



iH Q u^ = H = uj {b^b - a 2 )a 



The transformation of Hi creates a complicated expression H[ that connects states with different phonon 
numbers. This makes it difficult to perform the iteration of the recursion relation. For Pq being the 
projection on the Hilbert space with N — phonons after the Lang-Firsov transformation was applied, 
Eq. © gives 

P (z - H'r'Po = [Po(z - H )P - P H[Pi(z - H'^PiH'iPo]' 1 , 
and by using the approximation (z — H')^ 1 ps (z — Hq)^ 1 one obtains 

Po{z - HT'Po w [PoO - H^Po - P a H[Pi{z - H^PiH^a]' 1 . 
After a lengthy but straightforward calculation this becomes 

-Yi-z/LUo-a 2 ,-^ 2 )] 



t 2 e-° 
u V (z/ojq + a) 

with the incomplete Gamma function^ 

7*0, y) = 

m>0 



CO 



(7) 



i (-y) m 

m! a + m 



The renormalization factor e~ Q ' of the tunneling rate t is a well-known effect of the phonons, originally 
established in polaron physics^, and also observed in the strong-coupling regime of the Hubbard- Holstein 
model 1 -. Results of the iteration are shown in Fig. 1 for the spectral density. It should be noticed that the 
expression for Pq(z — H)~ 1 Po with the projection Pq before the Lang-Firsov transformation was applied 
is more complicated. 
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B. Weak Electron-Phonon Interaction 



If the tunneling energy is large in comparison with the electron-phonon interaction a different separation 
of the Hamiltonian H — Hq + Hi is needed: 

Hq = ta\ + Wob'bao, Hi = g<r a (b^ + b). 

Eq. (J3J can be truncated for N phonons. This leads to the equations 



G 



2,7 



za Q - tai - j(T Q - g 2 a 3 bG 2 j+ 2 b^a 3 



-l 

N-j 



with terminating condition 

G 2N = {{z - N)a - taiY 1 = 



I 



z-N 
(z -N) 2 -t 2 V t 

Some results of the iteration for the spectral density are shown in Fig. 2. 



z-N 



C. Two sites with two electrons 

The Hamiltonian reads in the basis of Eq. JSJ 



H e = 



( -t -i s 

t t 

-t t U 

\-t t u 



H ph = uj (b\bi + b\b 2 ) 



/I 











1 











1 


Vo 









and 



H ep h — 9 



/-At -A CP 

At + A 



V o ooo, 



7(At + A)5, 



(8) 



(9) 



where A — bi — b 2 - Thus A (A T ) lowers (raises) the number of phonons by one. 

Using Hq = H e + Hph and Hi = H ep h (this is weak-coupling case) the recursion relation of Eq. J1J 
reads 



Ga 



Pj{H e + H ph )Pj - giPjASGj+iSA^Pj 



with the terminating condition for a maximum of N phonons: 



G 



N 



z-P N (H e + H ph )P t 



N 



N 



(10) 



(11) 



This A(N + 1) x 4(7V + 1) matrix is diagonal in terms of the phonon states. There are N + 1 different 
phonon states, since there can be k (= 0, l...,N) phonons at the first site and N — k phonons at the 
second site. 

The action of the phonon operators PjA and A^Pj on the 4(j + 2) x A(j + 2) phonon-diagonal matrix 
Gj+i creates a 4(j + 1) x 4(j + 1) phonon-diagonal matrix matrix with 4x4 matrices D(k,j + 1 — k) = 
G J+ i(k,j + l-k\k,j + l-k): 



P 3 A 



/D(0,j + 1) 








V 







A T P.i 



D(j + 1,0) / 



G 



/D(0,j) ... \ 



V ... D(j,0)J 

This can be used to perform the iterations according to Eqs. (|10|l and 111(1 . The results for the spectral 
density of the electronic states | | d) and | |J,, 0) with a maximum of -/V = 3 phonons are shown in Fig. 
3. 



VI. DISCUSSION 



The recursive evaluation of the projected Green's function of Sects. 5.2 and 5.3 can be easily performed 
with an algebraic manipulation program. To demonstrate the qualitative tendencies the case for small 
numbers of phonons N = 1,2, 3 is discussed in the following. It should be noticed that much higher 
numbers can be studied with little effort. They show the same tendencies as those of small N. 



A. Single-site model 

In the strong-coupling regime the tunneling between the orbitals can be completely suppressed (i.e., 
t = 0). Then the projected Green's function has only one pole, namely z = 0. For any t > 0, however, 
there are infinitely many poles due to the incomplete Gamma function in Eq. (|7|). But not all poles 
contribute with the same weight, as it is shown in Fig. 1. There are two effects: (i) Phonons lower the 
groundstate energy and (ii) the weight of the excitations decreases rapidly with increasing energy. 

Eigenvalues in the weak-coupling regime in the absence of phonons (i.e. N = 0, where Hq = to\ + 
ujob^bao) are E — ±t. This is a level splitting caused by the tunneling between the two orbitals. Already 
a single phonon lowers the groundstate and creates new excited states, as shown in the second plot of 
Fig. 2. Additional phonons do not affect the groundstate but shift excited states and create new ones. In 
Fig. 2 this is plotted for phonon numbers up to N = 3. The spectral weights of the excited states are 
also affected by the increasing number of phonons. All these effects are related to fact that the elements 
of the projected Green's function are rational functions 

Pn{z)_ 

Qn{zV [ ' 

where N is the maximal number of phonons taken into account in the virtual processes. For t = g = 1 
and N = 0, 1, 2,3 the Green's function ((t \{z - t) has been calculated with MAPLE as 

N = 0: 



1 



N = 1: 



-4 z + z A + 7 z 2 + z 5 - 4 z 4 - 1 
-5 + 4 z 2 + 12 z + z A - 14 z 3 + 4 z 5 - z 6 

4 + 3 z - 46 z 2 + 33 z 4 + 24 z 3 + 10 z 6 ~ 33 z 5 - z 7 
~ -11 - 64 z + 46 z 2 + 94 z 3 + 31 z 6 - 14 z 5 - 82 z 4 - 10 z 7 + z s 

8 - 218 z + 701 z 2 - 558 z 3 - 449 z 6 + 759 z 5 - 361 z 4 + 128 z 7 + z 9 - 18 z 8 
~ -180 + 872 z - 716 z 2 - 1156 z 3 - 519 z 6 - 440 z 5 + 1809 z 4 + 414 z 7 - z 10 + 18 z 9 - 126 z 8 

Apparently the order of the polynomial P/v(z) is one less than the order of Qn(z) and increases by 2 for 
N > 1: 

P N (z) = z 2 ( N+ V +1 + (z 2 < Ar + 1 )), Q N (z) = z 2 ^ + oiz^+V- 1 ) 



N = 2: 



N = 3: 
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B. Two-site model 



In the case of two electrons on two sites there is tunneling between the two sites. Without phonons 
the Hamiltonian Hq = H e in Eq. JSJ has eigenvalues 



E = 0, U, C/± ^ TW 

Without Hubbard repulsion U there are three different eigenvalues, where the degeneracy of E = 
represents the states 

E = 0: \ U) + | U), E = U : \ ||,0) + |0,J.T). 

With Hubbard repulsion U > there are four non-degenerate eigenvalues. The most obvious effect in 
this two-site model is the lowering of the groundstate energy already by a single phonon and the creation 
of excitations on energies above the groundstate energy. It is interesting to notice that more excitations 
contribute to the singly-occupied Green's function (J.,T \(z — H)^ 1 ] than to the doubly-occupied 
Green's function (|t,0|(z — 4T,0). Moreover, there is always a state with maximal weight at 

E = U for the state | J.T,0). 

The Green's function with maximally N phonons is again a rational function of the form 1|12[1 . For 
t = g = 1 and U = the expression 

U,T|(z-i/)- 1 |4,T) = ^4 

WN\Z) 



has been calculated with MAPLE as 
N = 0: 



N = 1: 



N = 2: 



P (z) = z 2 -2, Q ( 2 ) = z(z 2 -4) 

P t ( z ) = ( z 5 - 3z 4 - 7z 3 + I5z 2 + 8z- 10)(z 2 - 2z - 5) 
Qi(z) = z 8 - 5z 7 - 10z 6 + 62z 5 + 33z 4 - 221z 3 - 44z 2 + 176z - 40 

P 2 (z) = z 10 - 13z 9 + 49z 8 + llz 7 - 378z 6 + 386z 5 + 720z 4 - 852z 3 - 536z 2 + 408z + 144 



Q 2 (z) = z 11 - 13z 10 + 45z 9 + 61z 8 - 562z 7 + 396z 6 + 1788z 5 - 2092z 4 - 1368z 3 + 1632z 2 + 224z - 192 
N = 3: 

P 3 (z) = z 13 -24z 12 + 226z 11 -1000z 10 + 1578z 9 + 3192z 8 -15756z 7 +15484z 6 +16149z 5 -35976z 4 + 8226z 3 

+15076z 2 - 7944z + 672 

Q 3 (z) = z(z 13 -24z 12 + 222z n -906z 10 + 708z 9 + 7046z 8 -22638z 7 +9330z 6 + 61467z 5 -98222z 4 -|-18512z 3 

+56152z 2 - 38816z + 7168) 

Apparently, the order of the polynomial Pat(z) is one less than the order of Qn(z) and increases by 3 for 
A - !: 

Pv(z) = z 3 < Ar + 1 ) +1 + o(z 3 ( w+1 )), Q N (z) = ^( JV +D+ 2 + (z 3 < Ar+1 ) +1 ). 



VII. CONCLUSIONS 



A a recursive method has been developed to study the spectral properties of small Jahn- Teller systems. 
It is based on a decomposition of the infinite-dimensional Hilbert space, spanned by a few electronic 
and an unlimited number of phononic states. Two cases were considered, one for an infinite number of 
phonons (the strong-coupling case) and one in which the number of phonons is increased by one in each 
iteration step of a recursive equation (the weak-coupling case). In both cases the iteration of the recursion 
relation leads to a continued-fraction representation of a projected Green's function. The matrix elements 
of the resulting Green's function are related to the incomplete Gamma function in the strong-coupling 
case and to rational functions in the weak-coupling case, respectively. 

An advantage of this method is that it approaches the exact solution systematically by standard 
functions. It is an alternative to perturbative approaches, based on a power series of a model parameter 
(e.g. the tunneling rate t or the electron-phonon coupling constant g). 
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FIG. 1: Strong coupling approach: The spectral density for z = E + 0.05i and g = 1, N = 1, and t = 0, 0.5, 1. 
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FIG. 2: The spectral density of the projected Greens's function for the single-site model with t = g = 1 and 
N = 0, 1: — Im((1 \{z — H)~ 1 \ f) (upper curves) and 7m({| |(jz — 1} (lower curves). 
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FIG. 3: as previous Fig. but for N — 2,3. 
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FIG. 4: The spectral density of the projected Greens's function for the two-site model with t = g = 1, U = 0, 
andiV = 0, 1: -7m((4,T |,T) (upper curves) and Jm((jT, 0|(z - #) _1 | |T, 0) (lower curves). 
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FIG. 5: as previous Fig. but for N — 2,3. 
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FIG. 6: as previous Fig. but for U = 2, N = 0, 1. 
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FIG. 7: as previous Fig. but for N — 2,3. 



